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<]} [ Abstract 

' The low energy structure of a theory containing light and heavy 

particle species which are separated by a mass gap can adequately be 
^ , described by an effective theory which contains only the light particles. 

^ | In this work we present a thorough analysis of the effective field theory 

of the linear O(N) sigma model in the spontaneously broken phase. In 
particular, we present a detailed discussion of two techniques, a short- 
distance expansion and a method based on loop- integrals, which can 
be used to explicitly evaluate the functional relationships between the 
low energy constants of the effective theory and the parameters of the 
underlying theory. We furthermore provide a detailed analysis of the 
matching relation between the linear sigma model and its effective 
theory, in order to clarify some discrepancies which can be found in 
the literature. 
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1 Introduction 



The effective field theory of the linear sigma model has recently received 
considerable attention in the analysis of the symmetry breaking sector of the 
electroweak interaction. So far the agreement between experimental data and 
the standard model predictions is almost perfect. Yet, the precise nature of 
the spontaneous breakdown of the electroweak symmetry remains unknown 
and various scenarios are discussed in the literature, including such different 
theories as a heavy Higgs |I|, |2|, |3|, || or technicolor models H . One way out 
of this ignorance is to replace the Higgs sector of the standard model by a 
model independent parametrization which provides a unified description of 
the low energy physics in the symmetry breaking sector, i.e. up to energies 
of the order of 1 TeV. If all the unknown particles of the underlying theory, 
e.g. a Higgs or a Technirho, are heavy, physics at low energies is dominated 
by the Goldstone bosons arising from spontaneous symmetry breaking, and 
all the other known particles of the standard model. In this case a con- 
venient description is provided by replacing the linear sigma model in the 
standard model Lagrangian by an effective Lagrangian which contains only 
the Goldstone bosons [p], |(|. The unknown physics of the symmetry breaking 
sector is then hidden in the low-energy constants which occur in the effective 
Lagrangian. Thus, if the relations between the constants in the effective La- 
grangian and the parameters of the underlying theory are known, low energy 
physics can provide interesting information about the precise nature of the 
spontaneous symmetry breaking which may occur at a comparatively high 
energy scale, such as in technicolor models. 

In this work, we present a thorough analysis of the effective field theory de- 
scription of the linear O(N) sigma model in the spontaneously broken phase. 
Our goal is twofold. On the one hand, we provide a detailed discussion of the 
matching relation between the linear sigma model and its effective theory. A 
formal definition of the effective theory as an adequate representation of the 
full theory in the low energy region will certainly require that corresponding 
Green's functions in both theories have the same low energy structure. It 
is known |7], |J that literally integrating out only the heavy degrees of free- 
dom in the full theory generally does not yield an effective Lagrangian which 
meets this requirement. Our discussion of this point is intended to clarify 
the discrepancies which can be found in the literature on the effective 
theory of the linear sigma model. On the other hand, we want to go beyond 
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a formal definition and provide a detailed discussion of two techniques which 
can be used to explicitly evaluate the effective Lagrangian. The first one em- 
ploys a short distance expansion in configuration-space, while the second one 
is based on the evaluation of loop-integrals in momentum-space. For the case 
we are considering here the result for the effective Lagrangian can already 
be found in the literature ||. The derivation given there, however, is rather 
sketchy and does not discuss the steps involved in any detail. The present 
article contains a thorough account of the problem. Furthermore it focuses 
on the discussion of the general techniques which can be used to evaluate 
the effective Lagrangian for a given underlying theory, provided the theory 
shows a mass gap between heavy and light particles and the coupling is weak 
in the low energy region. 

The outline of this work is as follows: In the next section we present 
a brief review of the linear sigma model. In section 3 the general low en- 
ergy structure of the effective Lagrangian is discussed, exploiting only the 
symmetry properties of the full theory. In the following section we provide 
a thorough analysis of the matching relation between both theories which 
includes an exact definition of the effective Lagrangian. In the next three 
sections the detailed description of two methods to calculate the low-energy 
constants at order p 4 can be found. In section 8 we discuss renormalization 
and express the bare quantities in terms of physical parameters. Finally we 
summarize this work in section 9. 

2 The Linear 0(N) Sigma Model 

In this section we employ the technique used in ref. from which the abun- 
dant literature on the linear sigma model may be traced. 

In the absence of external fields the Lagrangian of the O(N) symmetric 
linear sigma model is given by 

\d»fd^ + \m 2 f<t>-\(f<tf , (2.1) 

where (p A is an iV-component field. For m 2 > the classical potential has 
its minimum at a nonzero value T = m 2 /g and the O(N) symmetry is 
spontaneously broken to 0(N — 1). Accordingly, the iV-component field <p A 
describes 1 massive (pseudo) scalar and N — 1 massless Goldstone bosons. 
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In order to obtain the generating functional for the Green's functions we 
couple the field <fi A to a set of external fields with spin and spin 1: 

C a = Uj^ T V^ + \m 2 f<\> - | (ftf + f T <f) + htxF fll/ F^. (2.2) 

The spin 1 field is defined in terms of the generators T a of the Lie algebra 
o(N) and the N(N - l)/2 gauge fields F£, 

F AB = {Ta) AB F a ^ (2 3) 

which are coupled covariantly to the scalar field, 

= d^ A + F AB <p B . (2.4) 
The generating functional W a [F^, f] is defined by a path integral 

e«W./l = [dnWe'f* 1 **- . (2.5) 

Derivatives of this functional with respect to the spin field f A generate 
Green's functions of the scalar fields (p A while derivatives with respect to the 
spin 1 field F^ generate Green's functions of the currents, 

j; = (d^) T T> . (2.6) 

The last term in the Lagrangian is the trace over the square of the field 
strength F^f associated with the nonabelian gauge field F AB , 

F^u = [V M , V,] . (2.7) 

In the absence of the spin 1 field F^ the generating functional W CT [0, /] can 
be rendered finite by a proper renormalization of the mass m, the coupling 
constant g, and the scalar fields (p A and f A . Green's functions of the currents, 
on the other hand, are more singular at short distances than Green's functions 
of the fields. The time ordering of current operators gives rise to ambiguities 
which do not occur for the fields themselves. The corresponding Green's 
functions are unique only up to contact terms. This ambiguity is reflected 
by the presence of the last term in the Lagrangian (|2.2j ), which enters Green's 
functions of the currents through contact terms. The corresponding constant 
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h can be renormalized such that the full generating functional W^F^, f] 
remains finite after the regulator is removed. 

The generating functional defined in eq. (|2.5|) has two interesting proper- 



ties. First, note that both the generating functional W a and the Lagrangian 
C a are gauge invariant under local O(N) transformations of the form 

<t> - V<f>, 

f - Vf , (2.8) 
Ffj, -> VFpV T + VdpV T . 

This statement is equivalent to the Ward identities of the linear sigma model. 
Second, this description also allows for the presence of an explicitly symmetry 
breaking term of the form c T which provides the Goldstone bosons with a 
mass proportional to \J\c\. Keeping this constant small enough, one still 
has a mass gap between the Goldstone bosons and the heavy particle. To 
obtain Green's functions in this case one merely has to expand the generating 
functional around the nonzero value f — c instead of / = 0. Except for some 
minor changes our analysis remains valid for arbitrary values of the constant 
c, as long as a mass gap is maintained between the Goldstone bosons and 
the heavy particle. 

Let us pause for a moment and consider the particular case N = A which 
describes the symmetry breaking sector in the standard model of the elec- 
troweak interaction. The group 0(4) is semi-simple, 0(4) « SU{2) x SU(2), 
and one can express the nonabelian gauge field F^ in terms of external vector 
and axial-vector fields: 

K = < (2.9) 
pjd = _ e ijk v k_ (21Q) 

Using the method of steepest decent, the one- loop approximation to the 
generating functional is given by 

W a [F„ f}= J d d x£ CT (0 o , d^o, F v , f) + ^lndetD , (2.11) 

where D is the differential operator 

(y, Dy) = J d d xy T (V^V + a) y (2.12) 
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with 



~AB 



(-m 2 + gfy Q ) 5 AB + 2g^4>o • (2-13) 
The field <p A is a solution of the classical equations of motion 

V M V> = / + m 2 o - g($<f> )<j> • (2.14) 



The one-loop approximation given in eq. (|2.11 ) represents the first two 



terms in the expansion of the generating functional in powers of h. This 
corresponds to expansions of the Green's functions in powers of the coupling 
constant g. In the next step we will analyze the low energy structure of this 
approximation in the range of small energies and momenta, i.e., below the 
(physical) mass M of the heavy particle. This will in addition introduce an 
expansion in powers of momenta p 2 /M 2 . Since we want to describe the low 
energy structure of the linear sigma model in terms of an effective Lagrangian 
which contains only the Goldstone bosons, the following parametrization of 
the classical solution (p A will turn out to be useful: 

^ = —RU A , U T U = 1 . (2.15) 

The massless modes are now represented by the iV-component vector U A , 
confined to the sphere SV-i while the massive mode is described by the radial 



variable R. It is well known ||10|| , that the coordinates of the coset space 
related to the spontaneous symmetry breaking pattern, i.e. the full group 
factored by the unbroken subgroup, provide a convenient representation for 
the Goldstone boson fields in the effective Lagrangian. In our case the coset 
space is 0(N)/0(N — 1) which is isomorphic to the sphere Sjv-i- Thus, the 
definition ( |2.15D turns out to be quite obvious. 

In terms of the new field the minimum of the classical potential is located 



at R = 1 and the equations of motion Q2.14Q can be rewritten in the form 



DR + R(U T V,yU) = XqU + m 2 R(l- R 2 ) (2.16) 
R(V^ l U-U(U T V^U)) = xo - U(xo U) - Id^RVU , (2.17) 

where 

Xo = ^-f- (2.18) 
m 
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The radial variable R describes a massive fluctuation around R = 1 and 
its value at a particular space-time point depends only on the behaviour of 
the external field in a small neighbourhood of this point. Thus, for slowly 
varying external fields the behaviour of the massive mode is under control 
and the equation of motion ( |2.16| ) can be solved algebraically. The result is 
an expansion in powers of derivatives of the external fields: 

R = 1 + 5R l + 5R 2 + . . . , 5R n = 0(p 2n ) . (2.19) 

Note that the field / counts as a quantity of order p 2 while the gauge field F M 
occurring in the covariant derivative is of order p. U is a quantity of order 1. 
The first two nontrivial terms of this series are readily evaluated to be 

SRi = i^[xoU + V,U T WU] (2.20) 
2 1 2m 

The tree-level contribution to the generating functional is given by 

i~ 
9 



SR 2 = --8R\ + -^—^SR^V^VU) + total derivative . (2.21) 



\— J d d x (R(xoU) + \™ 2 R 4 ) +hj d d xtrF^F^ . (2.22) 



Thus, with the help of expansion (p. 19 ) one is able to determine the low 



energy behaviour of the tree-level contribution to the generating functional 
to any required order in p 2 . 

Finally, by changing to a new basis, we separate the massless fluctuations 
around the classical solution 0^, which are bound to the sphere U T U = 1, 
from the massive fluctuation along U A : 

N-l 

4> A ~ 0o = iU A + J2 <V . (2-23) 

i=l 

with 

efe j = 5 l3 , eJU = , (2.24) 
and with the completeness relation 

N-l 

U A U B + ]T ef ef = 5 AB . (2.25) 
i=i 
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In order to describe the action of the differential operator D in the new basis 
of the mass eigenstates we introduce the following notation: 



d = 


□ + 2m 2 - + 3m 2 (i? 2 - 1) 


(2.26) 


ft = 




(2.27) 


D, = 




(2.28) 


r*i _ 
I* 




(2.29) 


5 = 




(2.30) 


M T £) = 


(HO 


(2.31) 


D = 




(2.32) 


a ij = 




(2.33) 



One finds 

(y, Dy) = (£, dO + (e, 5rj) + (5 T £, rj) + (rj, Dij) . (2.34) 

To separate massless and massive fluctuations we furthermore diagonalize 
this expression with the transformation £ = £ + d^Srj: 

(y, Dy) = (£, di) + ( V , (D - 5 T d- 1 5) V ) . (2.35) 

Thus, we obtain the following representation for the one-loop contribution 
to the generating functional 

lndetl) = IndetD + lndetrf + lndet(l - D~ 1 5 T d~ 1 5) . (2.36) 



3 The Effective Lagrangian 

In this section we will analyze the general low energy structure of the linear 
O(N) sigma model. We will set up the effective Lagrangian description for 
the generating functional exploiting only the symmetry properties of the lin- 
ear sigma model. Thus, this discussion is completely general and applies to 
any underlying theory which has the same symmetry breaking pattern and a 
mass gap between the Goldstone bosons and the heavier particle species. To 
begin the discussion, let us recall that the generating functional W^F^, f] 
as defined in the previous section is gauge invariant under local 0(N) trans- 
formations as given in eq. (|2.8j). Due to the presence of Goldstone bosons 
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generated by spontaneous symmetry breaking, the structure of the Green's 
functions is nontrivial even at low energies. Only the contributions from the 



massive modes admit a Taylor series expansion. It is known [II] that, if 
the underlying theory (the linear sigma model in our case) does not contain 
anomalies and is Lorentz invariant, the low energy structure of the generat- 
ing functional can be described in terms of an effective field theory with a 
symmetric effective Lagrangian. 

In our case, this effective Lagrangian depends on an iV-component field 
U A confined to the sphere U T U = 1 which describes the N — 1 Goldstone 
bosons and which is coupled to the external fields f A and F AB . The lat- 
ter is coupled to the O(N) vector U A by the covariant derivative given in 
eq. ( j2.4|) . The effective Lagrangian is the most general functional invari- 
ant under Lorentz transformations and the gauge transformations given in 
eq. O, 

C eff = C eff (U, V,U, V„V V U, /,...) . (3.1) 

The Lagrangian C e ff is a sum of terms with an increasing number of deriva- 
tives, corresponding to an expansion in powers of the momentum, 

C eff = C 2 + U + U + ■ ■ ■ (3.2) 

where Li is of order p % . (If explicitly symmetry breaking terms are present, 
this expansion includes powers of the Goldstone boson masses as well.) The 
generating functional is defined as the path integral 

e iw eJf [F,j] = J dfi pyJd d *c eJ f _ (3.3) 

In general, there are two different kinds of contributions to the generating 
functional. On the one hand, we have tree-level contributions, given by 
the integral / d d x£ e ff, which has to be evaluated at the stationary point, 
i.e. with the solutions of the equations of motion. One the other hand there 
are contributions from loops, which ensure unitarity. General power counting 
arguments show []12" |, that n-loop corrections are suppressed by powers p 2n 
as compared to the tree-level. Thus, one-loop corrections with vertices of 
£2 are of order p 4 while those with vertices of £4 are of order p 6 , as are 
two-loop corrections. The suppression of loops by powers of p 2 which allows 
the perturbative study of the systematic expansion in powers of momenta is 
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related to the fact that Goldstone bosons do not interact at threshold. Thus, 
at order p A the generating functional W e //[-F)i, /] is given by 

Weff [F„ f]=J d d x (£ 2 + £4) + ~lndet£> + 0(p 6 ) , (3.4) 
where the term lndet-D describes the one-loop corrections with vertices from 

u 

The leading term of the generating functional is of order p 2 and given by 
the tree-level contributions of £2 

W eff [F„ f] = J d d xC 2 + 0{p A ) . (3.5) 

Imposing Lorentz invariance and taking account of the identities U T 'V l j,U = 
and U T V ^V^U = — V At t/ T V /i ?7, the most general effective Lagrangian at 
order p 2 which is invariant under local O(N) transformations involves only 
two low energy constants and is given by 

ciV^V^f/ + c 2 f T U . (3.6) 

In view of eq. (|2.18| ) it will turn out to be convenient to introduce another 
field \ A proportional to f A and define this Lagrangian in the form 

£2 = F 2 (i V^lFWU + X T U) ■ (3.7) 

Accordingly, we will replace the dependence of the generating functional W e ff 
and of the effective Lagrangian C e ff on the field f A by the dependence on 
the field \ A - The action in eq. ( |3.5| ) is to be evaluated with the solution U A 
of the classical equations of motion 

V^f/ - UiU^^^U) = X - U( X T U) . (3.8) 

Note that this solution coincides with the solution U A of eq. ( |2.17|) only at 
leading order. The corrections U A — U A are due to the massive mode and 
admit an expansion in powers of p 2 (see below) . 

To obtain the one-loop contribution of order p 4 we again use the method 
of steepest descent and parametrize the fluctuations around the solution U A 
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of the equations of motion in a way similar to the one used in the previous 
section on the linear sigma model, 



U A -U A = Cef~lcCU A + ... . (3.9) 

The vectors {U A , e A } are orthonormal, and thus obey relations as given in 
eqs. ( |2.24| , |2.25| ). The presence of the second term in eq. ( |3.9| ) is necessary 



to ensure that both vectors, U A and U A , are of unit length. Introducing a 
notation similar to the definitions (|2.26| - |2~33|) , 



ft = U T W^ t (3.10) 
= ^ + f M (3.11) 
fj" = efV.e-- (3.12) 

-a = _ r jm + (y^T WD + x T JJ y l ^ (313) 

the differential operator D in eq. ( |3.4|) turns out to be 

D = D^D" + a . (3.14) 

Finally, we have to specify the general effective Lagrangian at order p 4 . 
Note that the most general effective Lagrangian at this order is given as a 
linear combination of a maximal set of independent O(N) invariant terms 
of order p 4 . On the one hand, redundant terms can be eliminated by using 
relations of the form 

J <l d xV u V^U T V»V u U = Jd d xV^V li U T VuV v U + V v U T F v ' t V l JJ, (3.15) 

which are readily verified using the definition of the field strength given in 
eq. (|2.7D. On the other hand, the Lagrangian £ 4 contributes only at the 
classical level. Hence, the equations of motion ( |3.8|) can also be used to 
eliminate further redundant terms. From eq. ( P-8| ) one can infer the following 
two identities 

(x T V^U) + (x T U)(y,U T V»U) + (x T U) 2 = X T X (3.16) 
{VpWtFVvVU) - (V,U T V»U) 2 + ( x T u) 2 = X T X- (3.17) 
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The general effective Lagrangian at order p 4 contains eight independent low 
energy constants and can be brought into the form 



U = di(V \JJ T V»U) 2 + d 2 {V l JJ T V u U){V^U T V v U) 
+ d 3 ( X T U) 2 + d 4 { x T U){V»U T V»U) 
+ d 5 (U T F^ u F^U) + d & (y u U T F^V, L U) 
+ d 7 ( X T x)+d8tr(F, u F^) . (3.18) 



4 Matching 

In section two we have determined the generating functional of the linear 
sigma model up to one loop. In the next section we have described the 
low energy structure of this generating functional in terms of an effective 
Lagrangian which contains only the N — 1 Goldstone bosons. Exploiting the 
symmetry properties of the underlying theory one can reduce the number 
of low energy constants that are involved in the effective Lagrangian to two 
at order p 2 and to eight at order p 4 . All these low energy constants are 
determined by the underlying theory. To express the low energy constants 
in terms of the parameters of the linear sigma model we require that both 
theories yield the same Green's functions in the low energy region: 

W a [F,J] = W eff [F^ X ] ■ (4-1) 

Recall that we have absorbed one of the low energy constants in section 3 
by introducing the external scalar field x A - Thus eq. (|4.1| ) determines the 
relation between the scalar fields x A an d f A as well. Using eqs. ( |2.11| , f2.3(i| ) 
and ( |3.4D we obtain at order p 4 : 

J d d x C a + ^lndetD + ^lndetci + -lndet(l - D^ 1 5 T d^ 1 5) 

= J d d x (A + £ 4 ) + '-\ndetD . (4.2) 

We will solve this equation for the low energy constants in a couple of 
steps. To begin with, let us recall that both theories, the linear sigma model 
and the effective theory, contain Goldstone bosons. Thus, both sides of 
eq. (|4.2|) should reproduce all singularities associated with the light particles 
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and the corresponding contributions should cancel. This is indeed the case. 
At order p A all of these singularities are described by the determinants of the 
two differential operators D and D. A glance at eqs. ( |2.17|j3.8| ) shows that 
there is only a small difference between these operators, related to the fact 
that the fields U A and U A satisfy slightly different equations of motion. We 
will, however, see below, that the corresponding difference between lndet-D 
and lndet-D is of higher order and can be neglected in our analysis. Thus, at 
order p 4 our matching relation simplifies to 

J d d x C a + ^lndetd + ^lndet(l - D~ 1 5 T d~ 1 5) = J d d x (£ 2 + £ 4 ) . (4.3) 

Note that the two determinants on the left hand side of this equation are still 
nonlocal objects, while the right hand side is a local expression. The operator 
d, however, is related to the heavy particle in the linear sigma model and 
its determinant admits an expansion in terms of local quantities involving an 
increasing number of derivatives, which correspond to an expansion in powers 
of momenta. The second determinant on the left hand side of eq. ( |4.3| ) is a 
more complicated matter. Since it involves the operator D of the massless 
modes it is not a purely local object. The analysis in the following reveals that 
its nonlocality in the low-energy region shows up at the order p 6 . Thus, at the 
order p 4 all terms in the low energy expansion of the determinants in eq. ( |4.3| ) 
are local. The comparison of corresponding coefficients on both sides of this 
equation yields the functional relationship between the low energy constants 
and the parameters of the linear sigma model at this order. 

If one were to determine the low energy constants up to the order p 6 , the 
generating functionals of the linear sigma model and of the effective theory 
would have to be calculated up to the two-loop level. In this case, additional 
terms would occur which ensure that the matching condition contains only 
local terms up to the order p 6 . This will be explained in greater detail below. 

Before we use eq. ( J4.3| ) to determine the low-energy constants up to the 
order p 4 , some further comments on our matching relation ( |4.1| ) are in order. 
The crucial point to notice here is that the expansion of the nonlocal terms 
in eq. (|4.3|) takes place after all integrations over the fields U and R have 
been performed. More explicitly, this amounts to a relation between the full 
path integrals, i.e. 

J duMe'I^tWort = J d/#K/ dd ^^ . (4.4) 
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where C e ff = £2 + ^4 at order p 4 . To get a better understanding of this 
point, let us define the quantity T e ff[U, F^, f] (as in ref. 0) by performing 
the path integral on the right hand side of eq. ( |4.4j ) only over the massive 
particle, 

e V.ttPM = J dfi[R yJ^C„(R,U,F,J) ) (4 5) 

where we have used the parametrization ( |2.15|) for the field <p A as well as the 
identity 

d/i[0] = dfi[U]dfx[R] . (4.6) 
Our matching relation can thus be rewritten in the form 

J d/ift/y/^"^'*) = J Muy T ^^ F ^ . (4.7) 

Note that the quantity T e ff[U, F M , /] itself is a nonlocal object. Since it 
involves only the propagation of the heavy particle, it certainly admits an 
expansion in terms of local quantities as well. Such an expansion would 
introduce another density L e ff, defined by 

d d xL eff = r e// . (4.8) 

However, in eq. ((4.7|), which defines the effective Lagrangian £ e //, one first 
integrates over the massless modes U and then expands the resulting nonlocal 
contributions, as explained in the sequel of eq. (|4.3|) . One may now pose 
the question of whether it is permissible to change the order of these two 
operations. In that case one would first expand the integrand on the right 
hand side of eq. ( |4.7j ), and then perform the integration over the massless 
modes. If that procedure were correct, one would obtain the result that L e // 
and C e ff are the same and our matching relation ( ^4.7| ) between integrals 
could be replaced be an equivalent relation between integrands, i.e. 

e iJd^xC eJf (U,F^ x ) = J dfJi [ R ] e i J & d xC„{R,U,F^f) _ ( 4 _g) 

At the classical level this relation is obviously correct, since then 

J d d xC eff = J d d xL eff = T eff = J d d xC a , (4.10) 
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where C a is to be evaluated with the classical solution R of the equations of 
motion fl2.16| ). However, if quantum corrections are taken into account, it is 
generally wrong. Using the method of steepest decent, we again expand the 
field R as 

R = R + z, (4.11) 
where Rq satisfies the equations of motion (|2.16|) . The measure is given by 



dfi[R] = AfHR^dR = ArY[R$-\l + 0(z))dz . (4.12) 

X X 

Furthermore, in the dimensional regularization scheme the factor Rq' 1 in 
the measure does not contribute to the integral. Thus, at the order we are 
considering here, the measure is given by 

d/i[R] = Arj[dz , (4.13) 



and one obtains 



J d d x C a + ^lndetd = J d d x L eff . (4.14) 



Now one can compare the results for L e ff and C 2 + £4 as given in eqs. ( |4.14| ) 



and ( |4.3|) . At order p they both receive contributions from loop integrals 
which include only propagators of massive particles, represented by the term 
lndetd In addition to that, however, £ 2 + £4 furthermore receives contri- 
butions from loop integrals which include propagators of both, massive and 
massless particles, as described by the last term on the left hand side of 
eq. (4.3). These contributions are missing in the representation ( [4.14= ) of 



L e ff. Thus, if quantum corrections are taken into account, one must re- 
quire equality between the full integrals, as in our matching condition (|4.1|) . 
In general this requirement cannot be replaced by a relation between inte- 
grands. In other words, the integration over the massless modes and the 
expansion of the nonlocal objects do not commute. The relevance of mixed 
loops containing propagators of light and heavy particles is also discussed in 
the framework of QED in ref. |§. 
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5 Short-Distance Expansion 



In this section we will discuss the evaluation of the low-energy constants 
from eq. (|4.3|) with the help of short- distance expansions. The first term on 
the left hand side of eq. ( |4.3|) , which describes all tree-level contributions, is 
readily evaluated with the help of eqs. Q2.19| , p20| , |2.21|) 



J d d xC a = hj d d xti(F flv F^) 



+ ( y 1 / <' ! - 



-v u u T wu + xlu 



(5.1) 



+ 



d d x 



V»U T V»U) 2 + 2(V M C/' i 'V^)(x'o U) + (xoU) 



.Ttt\2 



The second term on the left hand side of eq. (|4.3| ) describes all those 
one-loop corrections which involve only heavy particles. We expand this 
contribution in powers of the difference om — d — (□ + 2m 2 ) which counts 
as a quantity of order p 2 . At order p 4 we get 



^lndetd = ^lndetd + ^Tr[rf V M ] - ^Tr[rf 1 a M ] 2 , 
with do = □ + 2m 2 . The second term is given by 



G M (0) J d d xa M (x) 



(5.2) 



(5.3) 



where Gm{z) is the Feynman propagator of a scalar particle with mass M = 
\/2m. The third term is of the form 

-~ J d d xd d ya M (x)a M (x+y)G 2 M (y) = -- J d d xa 2 M (x) J d d yG 2 M (y)+0(p 6 ) ■ 

(5.4) 

Note that only the local behaviour of the slowly varying external fields 
is relevant since the massive propagator cuts off large distances. Using 
eqs. (|2.20| , |2~2T|J2.26|) we obtain the following relevant contribution 



-lndetfi 



2m 2 A (2m 2 ) 
2A (2m 2 ) - 



d d x 
1 

16tt 2 



2(\/^U T ^U)+3(xoU) 

d d x iy^u T v^u) 2 
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3A (2m 2 ) + — ) J d d x {V,U T V»U){xlU) 
(JV2m 2 ) + ^)J d d x(x^Uy , (5.5) 



where 



1 M 2 

A (M 2 ) = Ao + ^ln— (5.6) 
Ao = ^-(^-1(^ + ^(1) + !)). (5.7) 



This leaves us with the last term on the left hand side of eq. ( |4.3|) which is 
the hardest to evaluate. If we again expand the localized quantity {x\d~ x \y) 
in powers of om — d — (□ + 2m 2 ), only the first three terms contribute at 
order p 4 



^lndet(l - D- 1 5 T d~ l 5) = --Tr[5D- 1 5 T do 1 ] 



+ t -Ti[5D- 1 5 T d 1 a M d 1 ] - '-Tr [5 D- 1 5 T d 1 ] 2 , (5.8) 

where do = □ + 2m 2 . Furthermore, for the last two terms only the leading 
singularity of the kernel (x\5D~ 1 5 T \y) , determined by the d'Alembertian, is 
relevant at this order. In particular, the relevant contribution to the second 
term is given by 

-2z / d d xd d yd d zfT(x)U(x + y)a M (x + z)d^G (y)G M (y-z)G M (z) 



= -2i / d'xft (x)Ux)a M (x) / d d yd d zd»d u G (y)G M (y - z)G M {z) + 0(p b ), 

(5.9) 

where Gq(z) is the Feynman propagator of a massless particle. Again the 
massive propagator cuts off large distances. Lorentz invariance requires this 
contribution to be of the form 

~il ddx tf ftl( ™ I dd y G M(y) ■ ( 5 - 10 ) 

In the same way, we obtain the following result for the relevant contribution 
to the third term in eq. (|5.8j) 

-At J d d x{flf y ){f p f a ) J d d yd d zd d vd^G (y)G M (y-z)d p d°G (z-v)G M (v) 
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d(d + 2) 
Using the identity 
we get at order p A 



{> 1 d d x (uinutn + 2(/jr)(/jn) / dv&(v) • 



(5.11) 

flU = V,U T V U U , (5.12) 



% -TY[5D- 1 5 T d Q l a M d Q 1 } - ^Ti[5D- 1 5 T d 1 ] 2 = 

- ^A (2m 2 ) + T^p^) / d d x (V»U T VU)(V,U T V„U) 

- (^3A (2m 2 ) + / d d x (V^ T V M C/)(^C/) . (5.13) 

The first term on the right hand side of eq. ( |5.8| ) is more complicated, 
because the next-to-leading singularity of the kernel (x\5D~ 1 5 T \y) is relevant 
in this case as well. The computation of the corresponding contribution, 
which may be based on either the short distance expansion or the calculation 
of two-, three- and four-point functions, can be found in the next two sections. 
Both methods yield the following result at order p 4 

^Tr^LT^do 1 ] = -2m 2 ^A (2m 2 ) - JLi) J d d x V ,U T ^U 

+ ^A (2m 2 ) - JLjjj / d d x {WU T V»U){V »U T V V U) 
~ (A (2m 2 ) - JL±) J d d x ( X T U)(WU T V,U) 

11/' / Trr\2 



d d x IT) 



(4tt) 2 12 

^ 2 ^ + jk-^)! AdxUTF ^ u 
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+ 



3 Ao(2m } ~ (4^36 



d d £ Xo Xo 



(An) 2 12 



(5.14) 



The expressions in eqs. ( |5.1|j5.5| , |3.13| , |5.14D provide the complete result for 
the expansion of the one-loop approximation to the generating functional 
of the linear sigma model in powers of the momentum up to the order p 4 . 
Now one can read off the relations between the low energy constants and 
the parameters of the linear sigma model. The contributions of order p 2 
occurring in eqs. ( p.l| , [5.5| , [5.14| ), which have to be matched with the effective 
Lagrangian £2 given in eq. ( |3.7|) , can be brought into the form 



d d x 



provided that we set 



F 2 = — 



m 
9 



-(V»U T V»U)+x T U 



g 12A (2m 2 



(An) 



X 



Xo 



l + g 6A (2m 2 )- 



(Anf 



(5.15) 

(5.16) 
(5.17) 



1 17- , 2x 1 35 

A (2m 2 ) - - — - — 

Ag 3 ° V ; 4vr 2 36 



The relation between x an d Xo then allows one to determine the remaining 
low energy constants by matching the contributions of order p 4 occurring in 
eqs. (|3.18|) with the effective Lagrangian £4: 

di 

d 3 
di 
d 5 



;A (2m 2 



1 11 

(4tt) 2 18 



1 11 

(An) 2 ~6 



T 9 " loX »< 2m2 ' - (4^ 



1? / 11 

-A (2m 2 ) - - — - — 
6 V ; 4tt 2 72 
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(4tt) 2 12 

d 8 = h . (5.18) 

Finally, we have to show that it was indeed consistent to neglect all 
corrections which are due to the difference U A — U A . By comparing the 
equations of motion for U A and U A given in eqs. fl2.17|JT8|) one infers that 
U A — U A receives contributions of order p 2 , due to the fact that R 1, 
and also of order g, since X 7^ Xo- First of all we can use the relation 
( |5.17|) to express xo by X everywhere in the generating functional W a . The 
replacement of U A by U A can lead to corrections of order p A only in the 
following contribution to the generating functional: 

X[U,F, X ] = Jd d x^-(\/,U T ^U + 2 x T U 



J d d x {y^v^u + x T u 



(5.19) 



In all the remaining terms one can safely replace U A by U A , since the corre- 
sponding corrections to the generating functional are of the order g and p 6 
which is beyond the accuracy of our calculation. In order to show that 

E[U, F, x ] = Z[U, F, x ] + 0(g) + 0(p 6 ) , (5.20) 

we make an intermediate step to separate the expansions in p 2 and in g. Let 
V A be the stationary point of the functional £[V, F, X ], subject to the condi- 
tion V T V = 1. Since U A itself extremizes the functional £[77, m/^fg, X o), it 
can easily be seen that U A — V A = O(g) and, therefore, 

Z[U,F,x} = W,F,x]+0(g) , (5.21) 

since at the stationary point only corrections proportional to the square of 
U — V can enter. Similarly, because U A extremizes the action of the nonlinear 
sigma model (cf. eq. ( |3.7| )) we have U A — V A = 0(p 2 ), from which we get 

nU,F,x] = nV,F,x] + 0(p 6 ) . (5.22) 
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Putting eqs. ( |5.21| ) and ( |5.22j ) together we have proven the claim in eq. ( |5.20| ). 
The main point in the whole analysis is that, at the stationary point, the 
replacement of V by V + 5V can only lead to corrections which are of second 
order in 5V . 



6 Employing the Heat Kernel 

In Euclidean space the short distance properties of the differential operator 
D, 

D = -D fl D l , + a, D M = 9 M + r M , (6.1) 
are governed by the Laplacian □. In d-dimensional Euclidean space one has 

(x\e xn \y) = (Aix\)- d ' 2 e~ z2 l AX , (6.2) 

where z = x — y. In order to determine the properties of the operator D at 

(6.3) 



short distances, we define |13 

-AULA 



x\e '"\y) = (47rA)- d/2 e~ z /4A H(x\X\y) 
The heat kernel H(x\X\y) satisfies the differential equation 



' d_ l 

d\ + A 



+ —z^D^ - D^D^ + a H(x\X\y) = 



with the boundary condition 



H(x\0\x) = 1 . 



Inserting the Taylor expansion 

oo 

H(x\\\y) = '£\ n H n (x\y) 



(6.4) 



(6.5) 



(6.6) 



n=0 



into the differential equation (|6.4| ), one obtains the following recursive rela- 
tions for the heat coefficients H n : 



1 + n + z M Df, ) H n+ i 



z^ Dp H Q 




0. 



(6.7) 
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Note that the covariant derivative which appears on the left hand side of 
the heat coefficients acts to the right on the space coordinate x. In the fore- 
going equation we have shown this explicitly by the arrows. In the following 

we will furthermore use the covariant derivative D, 



(6.8) 



which will always appear on the right hand side of the heat coefficients and 

act to the left on the coordinate y. Applying the derivatives and 
repeatedly to the relations in eq. (|6.7|) one obtains the following results for 
the heat coefficients 



H {x\x) = 1 (6.9) 

H^xlx) = -a (6.10) 

H 2 (x\x) = ^Y, v Y, u + l -a 2 - ^[D^D^a]} (6.11) 

H D u ) = (6.12) 

= ^ , (6.13) 



x=y 



x=y 



x=y 



Ho D v 



x=y 



where 



(6.14) 

The kernel (x\SD~ l S T \y) occurring in the traces of eq. ( p.8| ) can be rewritten 
in the form 



(D,f,) T + 2/J DA (xlD-'ly) (2 D v f„ + (D v f u 



(x\5D~ l 5 T \y) - 

which in turn can be expressed in terms of the heat kernel, due to the relation 

^ dA 



(xlD-M 



~ z / 4A H(x\X\ 



/o (47rA) d / 2 

It is convenient to define yet another kernel h(x\X\y) by 



(6.15) 
slation 

(6.16) 



{x\8D~ x 5 



-z 2 /4A 



o (47rA) d / 2 



h(x\X\y) 



(6.17) 
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If we use the following representation for the massive propagator 

dp ^i_9^2 
k (4vrp) d / 2 



Gm{z) = / Tj^b, e~~ 2m p , (6.18) 



the first term on the right hand side of eq. (|5.8|) can be brought to the form 

-i^-n_i r dXd p f Ad„Ad„, „-4a+±)-im> P 



Tr[^-vv] = ^y o ^r^m I dxd y e " TU+;J Ph ^ x \y) ■ 

(6.19) 

The kernel h(x\X\y) is readily determined from eqs. (|6.3| , |6.15| , |6.16| , |6.17l) 
and turns out to be 

h(x\X\y) = 2/J [5 P H^j {DM + 2(D p f p ) T (h D^j f v 
+ 4/J {5 P H D v \ j v + (D p f p ) T H(D u f u ) + ^/Jtf^ 



2/J [HD u )f u - 2/J A, /„ + flH{D v f v 



- (A^) T #/J - ^tf ■ (6.20) 

Note that in this equation all factors appearing on the right hand side of 
the heat kernel depend on the space coordinate y while all factors appearing 
on the left hand side depend on x. The exponential factor in the integrand 
in eq. ( |6.19j ) cuts off contributions from large values of z. Since the kernel 
h(x\X\y) is a smooth function for slowly varying external fields, it admits a 
Taylor expansion in powers of z. The expansion of the last term in eq. ( |6.20| ). 
for example, can be written in the form 

fl{x)H{x\\\y)f v {y) 
= (f?Hf v ) | + z p ((D p f p ) T Hf u + fT(D p H)f L 



x=y 



+ ( {D p D a f p fHf v + 2{D p f,Y\D a H)f v + f^(D p D a H)f v 

Z \ / x=y 

+ ... (6.21) 



If we furthermore insert the Taylor expansion ( |6.6|) for the heat kernel and 
keep only terms which are relevant at order p 4 we obtain 



h(x\X\y) = (D p f p ) T H (D u f u ) + 4/J (d p H 5 v ) /„ 
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^ "T zj t i o fT zj t i Z P ZcT 

A 



(DpD^Y HoU + fa D p D a H f. 



A 



(D p f p ) T H (D u f u ) - (D p DJ u ) T H Q f p 



+2fi[D„ Ha D„ ) ./',, - 2,/; (£>,£>„ #")//- 



f l TT f ~M ^ fTTTf 



2A 2 



{D p D a f P Y H f v + ft D p D a H f 



(6.22) 



Note that we have already omitted all terms that are odd in z or contain one 
of the factors given in eq. (|6.12|) because they will drop out anyway. Now 
one can perform the Gaussian integration over z and use the relations (|67 
to ( |6.13| ) for the heat coefficients to obtain (up to partial integrations) 



Tr[5D~ 1 5 T d [ 



d d x 



d\dp 



o (a + P y/ 2 



-2m 2 p 



(6.23) 



2/; r^/ w - 2/; + iDufv) + 
- jy+~p~j {( D M T ( D M + 2 ( D M T ( D »f») + 2 /J^A - #Vm} 



+ 



{{D p f u ) T (D p f v ) + {D^fipM + (D p f p ) T (D v f u )^ 



To bring this result into the form of eq. ( |5.14| ) we make use of the relation 
given in eq. Q5.12 ) as well as of the following identity 

(DJ u ) T (D p f a ) = V^ u U T V p V a U - (^ T V M V,f/)(f/ T V p V CT f/) . (6.24) 

Note that some of the terms in this equation do not occur in the effective 
Lagrangian given in eq. fl3.18p . These terms can be expressed in terms of 
the complete set of independent invariants that occur in the Lagrangian. In 
order to eliminate redundant terms we proceed as discussed at the end of 
section 3. Using the definition of the field strength F pv given in eq. (|2.7|) 



one can eliminate the various dependent terms of the form V^V^t/ VpV<j£7 
in favour of only V P V P U T V V V yU and some other terms, that are already 
present in the effective Lagrangian. Since we have furthermore used the 
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equations of motion of the effective Lagrangian to reduce the number of 
independent terms in £4, we should do the same in this well. From 

the equation of motion fl2.17| ) of the linear sigma model we infer a relation 
similar to the one given in eq. (|3.17 ), 



y.u'v.uy-xUo + ixouy 







(6.25) 



which holds up to corrections of order p 6 which can be neglected. Thus, using 
the identities ( 6.24 . B~25| ) to eliminate all terms of the form V y U T V P V a U 
we obtain the result given in eq. ( |5.14j ). 



7 Loops 

In the preceding two sections we have evaluated all one-loop contributions in 
the matching relation Q4.3|) which determine the low energy constants of the 
effective Lagrangian at order p 4 by means of the short-distance expansion in 
configuration space. In this section we want to compare this method with 
the evaluation by means of loop-integrals in momentum space. We shall see 
that at order p 4 this involves tadpole graphs and two-, three-, and four-point 
functions. Since we want to present a general comparison between these 
two methods - short- distance expansion versus loop-integrals - we will not 
discuss all contributions in detail. 

First, we will discuss the evaluation of the second term on the left hand 
side of eq. (|4.3|) which can be expanded as shown in eq. ( T5.2|) . The second 
term on the right hand side of this equation is a simple tadpole graph of 
the massive propagator as shown explicitly in eq. ( |5.3|) . The third term in 
eq. ( |5.2|) involves a two-point function with two massive propagators. In 
momentum-space it is given by (cf. eq. ( |5.4Q ) 

i f d d Pl r d^p 1 

" 4 J (^^iW(-Pi) J {2l ,y {M 2 - P 2 )(M 2 - ( P + Pl ) 2 ) ■ {7A) 

The loop-integral is readily evaluated and is regular at pi = 0. Hence, it 
admits an expansion for small momenta, 

r d d p 1 
J (2vr) d (M 2 -p 2 ){M 2 - (p + pt) 2 ) 
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_,^ o(MV _| pln ^_j +1 jj (7.2) 
-i f2A„(A/ 2 ) + -L") + 0(p?) , (7.3) 



where 



4M 2 , s 

At order p 4 only the leading term of the expansion of this two-point function 
is relevant. Thus one obtains the same result as given in eq. (|5.4| ). This may 
be compared with the short- distance expansion used in eq. Q5.4| ), where the 
integrand <tm{x)cm{x + y) is expanded in y. In both cases only the leading 
term is relevant at order p 4 . 

To discuss the evaluation of the contributions from the third term on the 
left hand side of eq. Q4.3|) , it is convenient to introduce the following short- 
hand notation for the n-point functions. We define the general two-point 
function by 

d ( 1 [ d d p p w .. .p^ n 

B l t 1 ...n n {p 1 ;m 1 ,m 2 ) = - / . . d . 2 ^—— 2 — f— ^ . 7.5 

% J (2n) d [mf — p 2 )(m2 — (p + Pi> > 

The function B(pi, mi, m-z) denotes the integral where no momenta appear 
in the numerator, e.g. B(pi;M,M) is given explicitly in eq. (|7.2p . Similarly 
we denote the general three-point function by C^ 1 , ^ n (pi,p 2 ] mi, m 2 , m 3 ) and 
the general four-point function by -D Ml ... Mn (pi,P2,P3; mi, m 2 , m%, m 4 ). 

Whereas the kernel (x|c/ _1 |?/) is related to the propagator Gm of the 
massive particle, the kernel (xl-D" 1 )?/), which enters in eq. (|5.8| ), leads to 
propagators Go of the massless particles in the loops. We can formally rewrite 
the operator D in the form (cf. eq. ( |2.32| )) 

D = □ (1 + □- 1 (( J D^) + 2T^ + a)) . (7.6) 

Therefore we get 

(x\(D- 1 Y k \y)=S tk G (x-y) (7.7) 
- J d d zG (x - z) ((D^) ik (z) + 2Ff(z)d^ + a lk {z)) G (z - y) 

+ 4 j d d zd d u V^{z)Yi k {u) G {x - z)d^G {z - u)d u G (u -y) + ... . 
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In the last two traces on the right hand side of eq. fl5.8| ) only the first 



term in expansion (7/7) is relevant. Note that this term also determines the 
leading singularity of the kernel (x\5D~ 1 5 T \y) . The relevant contribution to 
the second term on the right hand side of eq. (|5.8| ) is given by the integral 
(cf. eq. (|D) 

-2 / ^ ^ d f T (Pi + P2)r(-Pi)*M(-P2)C,M,P2; 0, M, M) , (7.8) 

which involves two massive and one massless propagator. The loop-integral 
which contributes to the third term in eq. ( |5.8j ) involves two massive and two 
massless propagators and is of the form (cf. eq. (|5.11| )) 

D pupa ( Pl ,p 2 ,p 3] 0,M,0,M) . (7.9) 

In expressions ( |7.8j ) and ( |7.9| ) the external fields account already for four 
powers in the momenta so that only the leading terms in the expansions of 
the function C pu (pi,p 2 ; 0, M, M) and D pupa (pi,p 2 ,p 3 ; 0, M, 0, M), which are 
of order 1, are relevant. Using the expressions for the leading terms as given 
in the appendix, one obtains the same results as given in eq. (|5.13|) . 

The evaluation of the first term on the right hand side of eq. ( |5.8| ) by 
means of the short distance expansion was more complicated, because in this 
case the next-to-leading singularity of the kernel (x\5D~~ 1 5 T \y) is relevant 
as well. In the language of loop- integrals this translates to the fact, that 
all three terms in expansion ( |7.7| ) of the kernel need to be taken 

into account. Hence, at order p A the first term on the right hand side of 
eq. ( |5.8| ) receives contributions from two-, three-, and four-point functions. 
In particular, one obtains 

\ I (^7 {b T (piH- P i)B( Pl ; 0, M) + 2ib T (p 1 )f p (-p 1 )B»(p 1 ; 0, M) (7.10) 
- 2tf^p 1 )a(-p 1 )B^p 1 ;0,M)+Af^p 1 )f u (-p 1 )B^(p 1 ;0,M)} , 



d d pi d d p 2 
{2n) d {2n) d ' 



ft (Pi + P2) ((D p ro)(- Pl ) + a(- Pl )) f v (-p 2 ) 
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+ b T ( Pl + P2 )T IM (- Pl )f u (- P2 ) 
- /J(Pi + p 2 )r,(-pi)a(-p 2 )] C^( PuP2 ; 0, 0, M) 
2<(pi + P2)T v (- Pl )f p (-p 2 ) [C^ p ( Pl , P2 ; 0, 0, M) 

+ PiC^( Pl , P2 ; 0, 0, M) + p?C^( Pl ,p 2 ; 0, 0, M)l ) , (7.11) 



and 



d d pi d d P2 d d P 3 rp 



(2-kY (2-kY (pi + ^2 + P3)r,(-pi)r p (-p 2 )/ CT (-p 3 : 



where 



D^( Pl ,p 2 , P3 -0,0,0,M) , 



b T = {D,n T + 2fT» 



The leading singularity of the first term on the right hand side of eq. 
determined by the d'Alembertian, i.e. 

- 2iTr[/" d, i £ /^o 1 ] . 



(7.12) 

(7.13) 
is 



(7.14) 



In expression ([7.1 0|) it is described by the term involving the two-point func- 
tion B^vipx] 0, M). Accordingly, the first two terms of the low-energy expan- 
sion of this two-point function are relevant at order p 4 . In this expansion 
the term of order 1 describes the contribution from the leading singularity 
while the term of order P 2 contributes to the next-to-leading one. The con- 
tributions from the heat coefficient Hi, occurring in eqs. (|6.22 ), are readily 
identified as well. They are given in expression (|7. 1 1|) by the term involving 
the quantity a. 

In this way one recovers the same result as was obtained by using the 
heat-kernel method. However, the evaluation of the term — ^Tr (5 D~ 1 5 T dQ 1 ) 
by means of loop-integrals is painstaking, since due to the formal inversion 
of the operator D in eq. ( [7.7[ ) the covariant form of the final result is not 
obvious. 

Finally, we should mention that all three- and four-point functions with 
two or more massless propagators have infrared singularities in 4 dimensions. 
They would have to be taken care of, if we were to extract the contributions 
of order p 6 from the loop integrals, using the next terms in the low energy 
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expansion of these functions. For instance, the next-to-leading order contri- 
bution of C^ u (pi,p2]0,0, M) would enter at this order. The analysis shown 
in the appendix reveals that the low momentum expansion of this function 
is of the form 

CW(pi,p 2 ;0,0,M) = P^{p 1 ,p 2 )+Q^{p 1 ,p 2 )]n{^)+highex order, (7.15) 

where P^ u and Q^ u are polynomials of second degree in pi and p 2 . Only P^ u 
contains a term of order 1, which leads to a contribution of order p A through 
eq. (|7. 1 1|) . However, terms involving logarithms, like the one in eq. (|7.15|) , 



do not contribute to the effective Lagrangian. 

Note that we already encountered such terms in our calculation. They 
were hidden in the nonlocal determinant IndetD, which cancelled in the com- 
parison of the generating functionals of the linear sigma model and of the 
effective theory in eq. ( f4.2|) . A similar situation would occur if we were to 
determine the contributions at order p 6 . Note that in this case one has to 
calculate the generating functional in the linear sigma model as well as in the 
effective theory up to the two-loop level. Therefore additional terms would 
occur beyond lndet-D, which represents the one loop graphs only. In this 
way one gets nonlocal terms in the linear sigma model and in the effective 
theory from loops, where a Goldstone boson propagates over long distances, 
leading to the above mentioned infrared problems. These nonlocal terms 
would, however, cancel in the matching of the generating functionals and the 
remaining local terms would lead to the coefficients di at order p 4 and to 
analogous coefficients at order p 6 . 

The short distance expansion is only appropriate to get these local terms. 
If we use, for instance, a representation of the propagator by the 

heat-kernel as in eq. (|6.16|) , it only makes sense to expand the heat-kernel 
for small values of the parameter A, which corresponds to small distances 
x — y. Infrared singularities on the other hand show up at large A where the 
expansion does not make sense, especially if massless particles are present. 

8 Renormalization 

In order to obtain finite results for physical quantities one has to renormalize 
the bare constants g, m and h before the regulator can be removed (there is 
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no wave-function renormalization at the order considered here, such that the 
external field f A does not require renormalization). This will indeed render 
the linear O(N) sigma model finite, because it is a renormalizable theory. 
To analyze the ultraviolet behaviour of the determinant of the differential 
operator D given in eq. (|2.12 ) we again employ the heat kernel method. In 



Euclidean space-time the (i-dimensional determinant may be defined as 



lndetD = - / ^Tre~ XD (8.1) 
Jo A 

= -(ATx)- d/2 J™ &\\- l - d/2 J & d x \xH{x\\\x) , (8.2) 



where the H is the heat kernel of the operator D defined in eq. ( |6.3| ). In this 
representation the divergences which are related to the ultraviolet behaviour 
show up at the lower end of the integration over A. For d — 0, 2, 4, . . . the 
determinant has poles. In order to identify the residues of these poles, we 
split the integration over A into an integral from to Ao and a remainder. 
Using the Taylor series expansion (|6.6|) we obtain 



-lndetD = - / d d x trH (x\x) + - — -— / d d x tvHi(x\x) 
2 d J d — 2 47r J 

1 1 



+ / d d x tiH 2 (x\x) + . . . (8.3) 



With the result given in eq. ( 6.11| ), the the pole-term at d = 4 turns out to 
be given by 



d 



Note that the differential operator D is positive in the O(N) symmetric phase 
of the linear sigma model. In this case the integrand in eq. (|8.1|) is expo- 
nentially damped as A goes to infinity, which expresses the fact that there 
are no infrared divergences in this phase. Since the divergence structure of a 
renormalizable field theory is the same in both, the symmetric and the spon- 
taneously broken phase, the pole is removed by the following renormalization 
prescriptions: 

g = g r (l-2(N + 8)g r (X (2m 2 ) + 6g)) + O(gf) (8.5) 
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m 



m 2 r (l 



2{N + 2)^ r (A (2m 2 ) + 8m 2 )) + 0(g. 



h 



h r + — A (2m 2 ) + 0(g r 



(8.6) 
(8.7) 



These relations introduce the five finite but otherwise arbitrary constants 
m r ,g r ,6m 2 ,Sg and h r . In the case of a (massive) Yang Mills theory, where 
the gauge field is a dynamical field, h represents the inverse square of 
the gauge coupling constant. In the present context, where is an external 
field, loops involving gauge fields do not occur. The renormalization of h en- 
countered here only accounts for the renormalization of the gauge coupling 
constant which are due to loops of the field cf> A . The finite renormalization 
constants 5m 2 and 5g are in general eliminated by requiring two renormal- 
ization conditions. In addition to that, we will also replace the remaining 
two independent parameters m r and g r by two other quantities, the physical 
mass M of the heavy particle and the low energy constant F, whose physical 
meaning will be discussed at the end of this section. 

For the rest of this section we will use Minkowski space notation. The 
mass M is determined by the pole position of the connected two-point func- 
tion 



—% 



5f°(x)5f°(y) 



W* [F„ /]] 



=/=o 



oro°(*)/(i,))|o 



(8.8) 



Note that we have embedded the unbroken symmetry group O(N-l) in such 
a way that it does not act on the first component 0° of the ^-component 
field (f) . Thus, (ft describes the massive mode while the remaining N — 1 
components % = 1, . . . , iV — 1, represent the Goldstone bosons. In order 
to evaluate the pole position, we use eqs. ( [2.11D and Q2.36j ) as well as the 
solution R of the equations of motion (|2.16|) , 



R(x) 



1 



d d yG M (x - y)x°o(y) 



(8.9) 



3m 2 / d d yd d vd d wG M (x - y)G M (y - v)G M (y - w)x (v)x (w) + ... . 



One obtains the following result for the pole position: 

1 



M z - iMT 



2mt 



2(N + 2)g r 5m' 



+ (N — 1)(1 + m) 



7T 



(8.10) 
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where M and T are the mass and width of the heavy particle. In terms of 
the renormalized quantities m r and g r , the low-energy constant F turns out 
to be 



. 9r 



1 + 2 9r(^ ~{N + 2)5m 2 + (N + 8)5g) 



(8.11) 



Thus, one can eliminate the two renormalization constants 5m 2 and Sg and 
replace the two parameters m r and g r by the physical mass M and the low- 
energy constant F by requiring 



M 



2m 2 r , 



p2 _ T Jh L 



1.12) 



Thus, 5m 2 and 5g are given by 



(N + 2)Sm 2 = 



327T 2 



v 9 V 1 v^, 



(iV + 8)^ = (N + 2)5m 2 - 



32tt 2 



AT 



The low-energy constants as well as the field x turn out to be 



(8.13) 
(8.14) 



M 2 



X = 



F 2 64tt 2 



./ 



3N - 10 A F 2 
«i = t. A + 



1 /ll6 + 9iV-27V37r 



d 2 
d 3 



6 

2 A _ 
3 96tt 2 



2M 2 192tt 2 
1 /ll nl M 2 \ 



M 2 

- (3iV - 10) In— 
fi 2 



N -1 



n 3 
F 2 



2M 2 64tt 2 



— (v*™-*^ -(JV-l) 1^' 
3 /i 2 



rf 4 = (iV-2)Ao + ^-^(l5 + iV-3v / 37r-(iV-2)ln^ 



1 , M 2 ~ 
+ ln^- 
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«7 — n 

1927T 2 

1 1 M 2 

* = 12 A » + ft ' + 384^ V' (8 ' 15) 

Note that the low-energy constants are of the form 

di^SiXo + di, (8.16) 
where the di are finite and expressed in terms of physical quantities. It is 



furthermore important to note that the matching relation (4.1) between the 
linear sigma model and the effective theory determines both contributions 
to the low-energy constants, the finite parts di as well as the poles in d = 
4. These pole terms, ^A , occur in tree-diagrams of order p 4 and exactly 
cancel the corresponding poles in the contributions from loops with vertices 
from £2, as represented by the determinant in eq. (|3.4|) . Since the effective 
theory describes the low-energy structure of the linear sigma model, the 
renormalization of the latter must indeed render the former finite as well. 
Let us now compare the results for the low-energy constants given in 



eqs. (|5.18|) and ( |3.15|) with those of refs. @. [|, [J|. Our results in eq. (|5.18|) 



where the low-energy constants are expressed in terms of bare quantities, 
agree with those in eq. (B.ll) of ref. [[|. However, the renormalization con- 
dition used in that article differs from the one given in eq. ( |8.12j ) above, which 
has to be taken into account if the results of our eq. ( |8.15| ) are compared with 
eq. (B.12) of ref. §. 

The authors of ref. integrate only over the heavy scalar field in the 
linear sigma model and, thus, obtain the functional T e ff as defined in eq. (|4.5| ) 
above. The expansion of this functional in terms of local quantities defines 
another set of parameters, d\ , which can be inferred from eqs. (|5.1| ) and (|5.5|) . 
We obtain, for example, d\ = d\ = d\ = 0, which agrees with the results of 
ref. j|. The constants d 2 , d 5 and d 6 , on the other hand, receive contributions 
from mixed loops and are different from zero. The effective Lagrangian, 



2 The low-energy constants di are related to those in ref. [[)) by h — (I3 ~ di,hi 
di + c?7, /12 = d$, and li — di for i = 1, 2, 4, 5, 6, and to those in ref. ^ by 4Li = d\, AL2 



-2Lq = d fi , and L 



10 
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however, which describes the low-energy structure of the linear sigma model 
involves the low-energy constants di as given in eqs. ( |5.18|) and (|8.15|) and 
not the parameters dj . 

The analysis of the low-energy constants in ref. || includes the effects of 
dynamical gauge fields as well. In that article, the contributions from mixed 
loops are taken into account and, furthermore, it was found that loops con- 
taining gauge fields do not contribute to the low-energy constants at order p 4 . 
For a more detailed comparison between the low-energy constants obtained 
in refs. || and || the reader is referred to ref. ||. 

Let us conclude this section with two remarks about the low-energy con- 
stant F. At the order we are considering here it is identical with the constant 
F^ which describes the coupling of the conserved currents to the Goldstone 
bosons. Recall that the gauge fields F^ couple to the conserved currents J° 
given in eq. (|2.(j|). By choosing a suitable basis for the N—l Goldstone boson 
states 1 7r fc ^> Lorentz invariance implies that these currents have the following 
matrix elements 

(0|j;(0)|7r fc (p)) = ^ fc F^, (8.17) 

which defines the coupling constant F n . Note that only the N — l currents 
which correspond to broken generators of the O(N) symmetry couple to 
the Goldstone bosons. According to definition (|8.17| ), the residue of the 
Goldstone boson pole of the two-point function 



5F?(x)5F?(y) 



-W a [F u , f] = (0\T(JZ(x) Jl (2/))|0) co?m +contact terms 

F»=f=0 

(8.18) 

is given by F%. Note that the contact terms do not contribute to the pole 
position. Since F n is a low-energy constant, it may as well be determined from 
our representation of the generating functional W e ff. The explicit calculation 
yields 

F = F n . (8.19) 

At tree-level, the coupling constant F n is identical with the vacuum ex- 
pectation value v of the field 0°, which is given by 



O|0 u (a:)|O> = JL^W a [F„f\ 



conn 



(8.20) 

F„=f=0 
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Beyond tree-level, however, both quantities differ by corrections of order g 
and the relation is given by 




(8.21) 

This reflects the fact that the constant F n is not renormalized, because the 
normalization of the currents is fixed by the commutation relations. The 
vacuum expectation value v, on the other hand, depends on the normalization 
of the fields and, hence, on a convention. 

Finally we want to make a comment on the validity of our first-order 
perturbative calculation. The physics of the linear O(N) sigma model in 
the spontaneously broken phase can adequately be described by an effective 
theory, which contains only the Goldstone bosons, if energy and momenta are 
small compared to the mass of the heavy particle, i.e. if p 2 /M 2 << 1. The 
perturbative expansion in the coupling constant g, on the other hand, breaks 
down if this mass is too large, since g ~ M 2 according to eq. ( 8.12 ). As an 
example we consider the Standard Model values N = 4 and v = 246 GeV. 
In this case perturbation theory in the parameter g can no longer be applied 
if the mass is of the order of 1 TeV. 



9 Summary 

In this work we have analyzed the low-energy structure of the linear O(N) 
sigma model from an effective field theory point of view. In the spontaneously 
broken phase the spectrum of this theory consists of one massive mode and 
N — 1 Goldstone bosons. At low energies, physical phenomena are domi- 
nated by the light particles of a theory, while the presence of heavier species 
manifests itself rather indirectly. In the low energy region, the nontrivial 
structure of the linear sigma model is related to the singularities of Green's 
functions which are associated with the Goldstone bosons. The effects of the 
heavy particle, on the other hand, admit a Taylor expansion and show up 
as simple power-like contributions. Hence, the low energy structure of this 
theory can adequately be described by an effective theory which contains 
only the Goldstone bosons. In the Lagrangian which describes the effective 
theory, effects of the heavy particle manifest themselves as an infinite tower 
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of nonrenormalizable interactions which are accompanied by an infinite se- 
ries of dimensionful low energy (coupling) constants. The requirement that 
both, the effective and the full theory describe the same low energy physics in 
principle determines the relationships between all these low energy constants 
and the parameters of the underlying theory. In more mathematical terms, 
the equivalence between the effective theory and the full theory requires that 
corresponding Green's functions have the same low energy structure. This 
can be expressed in a compact way as 

W a [F,J] = W eff [F„ X ] , (9.1) 

where W^F^, f] and W e ff[F^ x] are the generating functional of the linear 
sigma model and of the effective theory, respectively. They depend on a set 
of external vector fields F^ and scalar fields x or which are coupled to 
the matter fields. Derivatives of these functionals with respect to the vector 
fields Fp generate Green's function of the O(N) currents, while derivatives 
with respect to the scalar fields x or / generate Green's functions of the 
matter fields. Equation (|9.1| ) should not be understood as an identity but 
rather as an asymptotic equality in the low energy region. It determines the 
low energy constants as well as the relation between the external scalar fields 
X and /. More explicitly, eq. ( |9.1| ) requires the following relation between 
full path integrals 

J dfi[U]e l I ddc ^ u ' F ^ = J d f j,[U]d f j,[R]e i S dd ^ R ' u ' F ^ , (9.2) 

where the effective Lagrangian C e ff(U, F^, x) depends on the external sources 
and on the fields U of the Goldstone bosons, while the Lagrangian of the 
linear sigma model C a (R, U, F^, f) depends on the field R of the heavy degree 
of freedom as well. At the classical level, this relation between integrals is 
equivalent to the following relation between integrands 

e ifd*£ eff (U,F„ X ) = J dfl[R]e iJ^(R,U,F„f) 3) 

This equivalence, however, is lost if quantum corrections are taken into ac- 
count. In other words, beyond tree- level, the correct matching condition is 
given by relation (|9.2| ) between the full integrals. In this sense, integrating 
out heavy particles should rather be understood as replacing the full theory 
by an effective one which describes the same low energy physics. 
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At low energies, Green's functions receive nonlocal contributions which 
are related to the propagation of light particles over large distances. Since 
both, the effective and the full theory have the same light particle content, all 
nonlocal contributions drop out of the matching relation. The remaining local 
terms involve the propagation of heavy particles over short distances. It has 
been shown in this work how the heat kernel technique can be used to evaluate 
the short- distance expansion of these terms and thus determine the functional 
relationships between the low energy constants and the parameters of the full 
theory. These relationships can, however, also be evaluated by means of loop 
integrals in momentum space. In order to present a detailed description of 
both methods we have explicitly evaluated the effective Lagrangian for the 
linear O(N) sigma model in the spontaneously broken phase up to the order 
p A . In contrast to the short distance expansion, the evaluation based on 
loop-integrals is not covariant during intermediate steps of the calculation. 
The momentum space picture is useful, however, to understand why eq. ( |9.3| ) 
does not define the correct effective Lagrangian beyond tree-level. The local 
terms in the matching condition receive quantum corrections from two types 
of Feynman diagrams, cf. eq. ( |4.3|) , those which involve propagators only of 
the heavy degrees of freedom and those which involve propagators of both, 
heavy and light particle species. The latter ones are obviously missed in 
eq. ( |9.3| ). Finally we have expressed the low energy constants in terms of the 
mass M of the heavy particle and the decay constant F. 
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Appendix 

In this appendix we briefly want to discuss how the low energy expansion of n- 
point functions can be constructed. As a particular example we shall consider 
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the two three-point functions C fiu (pi,p 2 ; 0, M, M) and C flu (pi,p 2 ',0,0, M). 
Note, that the latter one has infrared singularities while the former one is 
regular for small momenta. In Euclidean space-time one has 



C IM ,(pi,p 2 ;0,m 2 ,m 3 ,) 
d d p 



(A.l) 



(27i) d p 2 ((p + pi) 2 + m%)({p + pi +p 2 ) 2 + ml) 
A convenient Feynman parametrization is given by 



a 



(1 - (3) ((p + Pl f + m 2 2 )+[3 ((p + Pl + p 2 f + 



(A.2) 

+ (1 - a)p 2 = p 2 + 2pk + A , 



with 



k = a{pi+f3p 2 ) 
A = a 



1 - (5) {ml +pi)+[3 (ml + ( Pl + p 2 f) 
The integration over p yields 

C M „ = £ daa £ dp (k^KF(0) + i ^F(l)) , 



where 



F(0 



1 r(3 - z - 1) 



(47r) d / 2 (A - k 2 f- l - d / 2 ' 
In the following we will furthermore use dimensionless momenta 

Pi 



M 



In the case of C fll/ (pi,p 2 ; 0, M, M), m 2 = m 3 = M and we have 

where the function 

f l (a,[3,q 1 ,q 2 ) = (1 - + fifa + q 2 f - a(q 1 + (5q 2 ) 



(A.3) 
(A.4) 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



(A.9) 
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is of order q 2 . Thus, the three-point function C lxl/ (pi,p 2 ; 0, M, M) is regular 
for small momenta and the low energy expansion is obtained by expanding 
the integrand in powers of fi . At leading order we get 



C^( Pl ,p 2 ;0,M,M) = - -<5 M „ A (M 2 ) + 



1 1 



32tt 2 2 



In the case of C^(pi,p 2 ; 0, 0, M), m 2 = & m 3 = M and we get 

A-k 2 = aM 2 F(a,[3, qi ,q 2 ), 

where 

F(a, [3, q u q 2 ) = (3 + (1 - /3)g 2 + /3(gi + q 2 f - a(?i + /% 2 ) 2 . 



(A.10) 



(A.ll) 



(A.12) 



Since F(a, (3, qi, q 2 ) vanishes iff (3 = and a — 1, we encounter an endpoint 
singularity in this case. Therefore the three-point function C^(pi, p 2 ; 0, 0, M) 
will have a branch point at q\ = 0. In order to analyze the behaviour of this 
function at small momenta, we separate the term linear in (3 and 1 — a as 
follows 

F(a, (3, q u q 2 ) = + g 2 ) + (1 - a)g 2 ) (1 + f 2 (a, (3, Ql , q 2 )) , (A.13) 
where the quantity 



f 2 (a,P,q!,q 2 ) 



P 



(3(l + q 2 ) + (l-a)q 2 



(-(3aq 2 2 + 2(1 - a) qi q 2 ) (A.14) 



is of order q 2 in the region of integration. Note, that this function is well- 
behaved even in the limit (3 — > + and a — > 1~. Thus we can expand this 
integrand in powers of f 2 . In this way one obtains 



fdaa C &(3Wk v F{$) 
Jo Jo 



(A.15) 



16tt 2 



11 



18 



tftf - - -lng 2 + -qgfi + -(tf «£ + 



+ 



and 



f 1 daa f 1 d(3F(l 
Jo Jo 

= -A (M 2 ) + 



(A.16) 



32tt 2 



+ 
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For completeness sake we finally present all the terms in the low energy 
expansions of the n-point functions that we use in section 7 (in Minkowski 
metric): 

B(p i; M,M) = -2A (M 2 )-y^ (A.17) 
B( Pl ;0,M) = -2A (M 2 ) (A.18) 
B"(p i; 0,M) = (VM 2 ) + ^£) (A.19) 
M 2 „„ / f 2 11 



S^(p i; 0,M) = -— <T A (M 2 )- 



32tt 2 2 
32tt 2 6/ 

2 „ ,, f2 , . 1 2\ 



+^<r (a (m 2 ) + 



3^H Ao(M) + 32^3j ( A - 2 °) 
C^( Pl ,p 2 ;0,M,M) = 1^ (a (M 2 ) + (A.21) 



C^( Pl ,p 2 ;0,0,M) = -<t(a (M 2 ) - _-) (A.22) 
C^( Pl ,p 2 ;0,0,M) = (a (M 2 ) - JL±) 



V ( A » (M > + 32^) 



permutations of (/iz^p) (A. 23) 

1 

12 



(^ + 3^4) (A.24) 



12 w " " " " " 7 

( A »< M2 > " d^) ■ < A - 25 » 
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